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Abstract 

We compare the divergence of orbits and the reversibility error for discrete time dynamical 
systems. These two quantities are used to explore the behavior of the global error induced by 
round off in the computation of orbits. The similarity of results found for any system we have 
analysed suggests the use of the reversibility error, whose computation is straightforward since 
it does not require the knowledge of the exact orbit, as a dynamical indicator. 

The statistics of fluctuations induced by round off for an ensemble of initial conditions has 
been compared with the results obtained in the case of random perturbations. Significant 
differences are observed in the case of regular orbits due to the correlations of round off error, 
whereas the results obtained for the chaotic case are nearly the same. 

Both the reversibility error and the orbit divergence computed for the same number of 
iterations on the whole phase space provide an insight on the local dynamical properties with 
a detail comparable with other dynamical indicators based on variational methods such as the 
finite time maximum Lyapunov characteristic exponent, the mean exponential growth factor 
of nearby orbits and the smaller alignment index. For 2D symplectic maps the differentiation 
between regular and chaotic regions is well full-filled. For 4D symplectic maps the structure of 
the resonance web as well as the nearby weakly chaotic regions are accurately described. 



1 Introduction 



The discrete time dynamical systems are widely used because they present a rich structure (regular, 
chaotic, intermittent and other types of orbits) in one or two dimensions and because the numerical 
evaluation of the orbits is straightforward. A large number of dynamical tools known as indicators 
of stability allow us to improve our understanding of dynamical systems: Lyapunov Characteristic 



2010) as well as Return Times Statistics (Kac, 1934), (Gao, 1999) 



Rosenstein et al. 


1993 


), ( 


Skokos 


, (Hu et a/., 2004 


I , (Buric et al. 



2005]). In the recent past many other indicators have been introduced not only to address the same 
problem of quantifying the degree of chaoticity of an orbit but also to perform the task fast. The 
Smaller Alignment Index (SALI), widely described in Skokos et al. (2002) and Skokos et al. (2004), 
allows to discriminate regular from chaotic orbits. Similarly, the Generalized Alignment Index 
(GALI), introduced in Skokos et al. (2007), is a family of highly efficient algorithms. Besides, the 
Mean Exponential Growth factor of Nearby Orbits (MEGNO) discussed in Cincotta et al. (2003) 
and Gozdziewski et al. (2001) is a quantity that gives a fast identification of the chaoticity of the 
orbit while its average slope estimates the maximum LCE (mLCE), see Maffione et al. (2011) for a 
test of the MEGNO. Fidelity and correlation decay are also tools that can be successfully used to 



Vaienti et al. ( 


2007) 


and 


Turchetti et al. 


(2010b 


, (Robutel & Laskar 


2001). 



As it is already known, in any numerical computation of a given trajectory, there is a round off 
error, and it would be interesting to study its relation with the chaoticity of the orbit and if possible 
determine its effect with the help of some dynamical indicator. 

The error between an exact and a numerical orbit is due to the finite precision used to represent 
real numbers and to the arithmetics with round off. This is unavoidable because the length of 
the binary strings representing real numbers must not change after arithmetic operations. The 
shadowing lemma is often invoked to state the existence of a true orbit close to a numerical orbit for 



chaotic systems (Katok & Hasselblatt 1997), (Hammel et al. 1987), (Chow & Palmer 1992 ), (Chow 



& Palmer, 1991). However it does not provide information on error growth for a given numerical 
orbit and the case of regular numerical orbits is not covered by such a lemma. The global error 
between the exact and the numerical orbit is unknown because the first one is not computable. 
Nevertheless an estimate can be provided by replacing the exact orbit with another one having 
very high accuracy. If the map is invertible the reversibility error can be computed without any 
reference to the exact orbit. Both errors have a similar behavior, namely an average linear growth 
for regular orbits and an average exponential growth for chaotic orbits and consequently can be used 
as dynamical indicators. Due to the correlation between the single step errors, there is a substantial 
difference with respect to the case in which the exact system is perturbed with random uncorrelated 
noise. 

This difference has been analysed in Turchetti et al. (2010b) by using the fidelity which measures 
the deviation of the orbits of a given map and its perturbation by integrating over all the initial 
conditions with the appropriate measure. In the case of regular orbits with random perturbations 
the decay of fidelity is exponential whereas with round off errors the decay follows a power law. In 
the case of chaotic orbits the asymptotic limit is approached super-exponentially for both situations 
(random uncorrelated perturbations and round off). The symplectic maps of physical interest are 
generally provided by the composition of a linear map with another one whose generating function is 
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the identity plus a function of position or momentum only. In this case the inversion is immediately 
obtained in analytic form. As most of the symplectic maps used in the literature are of this kind, 
the condition of invertibility of the map is not too restrictive. In addition, the reversibility error 
can also be applied to time-reversible systems of differential equations. 

Here we present an analysis of the reversibility error and its comparison with the divergence of orbits 
due to round off and other dynamical indicators of stability such as SALI, MEGNO and a finite time 
numerical estimation of the mLCE. Moreover we talk about the similarities and differences between 
an irreversibility due to the single precision round off and one due to the application of uncorrelated 
random noise in an orbit iterated with double precision. For two dimensional area preserving maps 
the reversibility error for a fixed number of iterations detects the various regions of phase space 
with different stability properties quite effectively as well as other dynamical indicators. Besides, 
the reversibility error allows to study the structure of the resonance web of a four dimensional 
symplectic map. 



2 Round off error methods 



In a computational device a real number x can be represented by a floating-point number x* 
according to Goldberg (1991) and working with base 2, can be written as 



that 



p-i 

A;=0 



e—k 



where dQ.di - ■ ■ dp-i is called the significand and has p binary digits dk, whose value is or 1, and 
where the exponent is an integer that satisfies Cuiifi < e < 

^max- The IEEE 754 standard states 
that for single precision p = 24, Cmm = —126 and Cmax = 127 while for double precision p = 53, 
&min — —1022 and Cmax — 1023. 

Consequently, a floating point number x* differs from the real number it represents and the 
relative error r^, defined by x* = x{l +rp), satisfies \rp\ < e = 

' " 0-24 



' p, d efined by x* - 
and Knuth (1973). Therefore, 



according to IEEE 754 we have that e 



as analysed by Goldberg (1991 ) 
: 2-^^for 



and 



single and double precision which roughly corresponds to 7 and 16 decimal digits, respectively. The 
arithmetic operations such as sums or multiplications imply a round off, which propagates the error 
affecting each number. Round off algebraic procedures are hardware dependent as detailed in |Knuth 
(1973). Unlike the case of stochastic perturbations, the error strongly depends on x. Suppose we 



are given a map M{x) then the error with respect to the numerical map M^:{x) after the first step 
is defined by: 



6i = e^i = xl-xi = M^{x) - M{x). 

Analogously, we define the local error produced in the n''^ step by 6n 
M(x;_i) where ^„ = ^„(x;_i). 
The global error 



(2) 



Gn = M:{x)-M^{x) 



(3) 
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accumulates all the local errors and explicit expressions can be written at first order in e. In the 
example of a regular map we take the translation on the torus defined as: 

M{x) = x + cj modi, (4) 

so M*(x) = X* + u!* mod 1 and the global error, which includes also the error to the mod 1 
operation, becomes: 

n 

G„ = e^ 6 = e(< + 'i^n), (5) 
fc=i 

where is a time average defined as the limit of Gn/n for n — >■ oo and Wn is a bounded 
fluctuation. 

For the chaotic Bernoulli map 

M{x) = qx mod 1, (6) 

we have that M^{x) = qx* mod 1 and that the single step and global error satisfy, respectively, 
\Si\ < Ceq and |G„| < Ceg". 

In order to compute the global error we need the knowledge of the exact map, which is usually 
precluded. A practical way to overcome this difficulty is to replace the exact map with a map 
computed with an accuracy 2^^ where P p. For instance, as p = 24 corresponds to single 
precision one might choose P = 53 corresponding to double precision. If p = 53 then one might 
choose P = 100 and so forth. There are available libraries which allow to compute with any fixed 
number of bits or significant decimal digits. The computation of the reference orbit is expensive 
if high precision is used, but there is no other way to evaluate the global error. As a consequence 
the "exact" orbit is achievable for a definite number of iterations which depends on P and the 
nature of the map. Taking into account what we have just mentioned, in the forthcoming numerical 
experiments, we will use the divergence of orbits, defined by: 

A„ = M^(x)-MS(x), (7) 

where Ms and Ms stand for single and double precision iterations respectively. 

If the map is invertible there is another option to overcome the difficulty of not possessing the 
true map. We define the reversibility error as 

Rn = M-^oM2{x)-x (8) 

which is non zero since the numerical inverse M~^ of the map is not exactly the inverse of 
M* namely M~^ o M=k(x) / x. Obviously the reversibility error is much easier to compute than 
the divergence of orbits (if we know explicitly the inverse map) and the information it provides is 
basically the same as the latter. Both quantities give an average linear growth for a regular map 
together with an exponential growth for a chaotic map having positive Lyapounov exponents and 
strong mixing properties. W^hen computing Rn we will set M* = Ms in order to compare with A„. 
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3 Variational methods 



In the forthcoming sections we will compare the performance of the indicators presented above with 
three well known and widely accepted dynamical indicators that are based on the behavior of the 
solution of the variational equations of the system. These are the finite time mLCE, the cumulative 
moving time average of MEGNO, and SALI. 

Let us briefly state them for discrete time dynamical systems of the form: 



f(Xn), (9) 



where Xn is the state vector at time n and f is a vector valued function. 

The concomitant discrete time variational equations, also called tangent map dynamics, associ- 
ated to a given orbit {x„}„gN are the following: 

v„+i = Df(x„) • v„, (10) 

where Df(x) is the Jacobian matrix of the function f and v„ is a deviation vector at time n. 



Skokos (2010) makes a historical review of the definition of the LCEs and its connection with 
the divergence of nearby orbits. He also states the theorems that guarantee the existence of the 
spectrum of LCEs and, in particular, the current definition of the mLCE in terms of the solution 
of the variational equations: 

mLCE= lim -ln||^, (11) 

n-s>oo n 1 1 Vol I 

with ||.|| some norm. For a chaotic orbit the mLCE is positive and this implies an exponential 
divergence of nearby orbits. On the other hand, for regular orbits mLCE is zero. 

In order to have a numerically computable quantity we define the finite time mLCE at time n 

as 

1 II II 1 " II II 
mLCE(n) = -lnp4 = -yin#^, (12) 
ll^oll l|vfc-i|| 



so that equation Eq. (11) can be reformulated as 

mLCE = lim mLCE(n). (13) 



n—^oo 



Cincotta et al. (2003) defines a biparametric family of MEGNO indicators: 



Ymjin) = {m + l)n^ V A;™ In rf^, (14) 



where m and j are integer numbers. They made experiments with 1^2,0) a-nd Yx-x and concluded 
that the last one allows both a fast classification between chaotic and regular orbits, and a clear 
identification of stable and unstable periodic orbits. Due to this fact, we will use Y = Y\-\ 
throughout the rest of this article. In order to reduce the fast oscillations that the time evolution 



of the MEGNO presents, in Cincotta et al. (2003) they use a time average of this quantity, namely: 
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1 

fc=i 



(15) 



Theoretically, the asymptotic evolution of Y{n) for any dynamical regime can be put into a 
single expression: 

-, , mLCE 

y(n) « — - — n + c, (16) 
where c ~ and 2 for chaotic and regular motion, respectively. 



We also compare our results with the SALI, introduced by Skokos (2001), that measures the 



degree in which a pair of initially linearly independent deviation vectors tend to become aligned. 
The underlying principle is that, for a chaotic orbit, a deviation vector under the tangent map 
dynamics changes in order to become aligned with the instantaneous most unstable direction. In 
other words, for almost every pair of initial deviation vectors (vq, Uq), the more chaotic the orbit the 
faster that the angle between them will reduce to zero. In the case of regular orbits, the behavior 
depends strongly on the dimensionality of the map: in maps with dimension > 4 the deviation 
vectors generally remain unaligned so the SALI tends to a positive non-zero value, while in 2D 



maps the two deviation vectors tend to align with a time rate that follows a power law. See ( Skokos 



et al., 2004) and (Skokos et al. 2007) for numerical tests of the SALI and its generalisation, the 
GALI family. Denoting the Euclidean norm with ||.|| the SALI is defined as: 



SALI(n) = min 



+ 



Uv 



Ur 



(17) 



To control the exponential increase of the norm of the vectors and avoid overflow problem, 

i every time step, keeping their norm equal to 1. 
if and only if the two normalised vectors have the 



Skokos et al. (|2004|) have normalised them, at every time step, keeping their norm equal to 1. 

[0, V2] and SALI 



Defined like this SALI(n) G 
same direction, being equal or opposite 



4 Error behavior on simple maps 

From the linear and exponential bounds on the errors it follows that the numerical orbit remains 
close to the exact one for a number of iterations proportional to 1/e in the regular case and to 
ln(l/e) in the chaotic one. 

We consider two types of models where the error grows linearly and exponentially respectively. 
The first one is the translation on the torus defined by Eq. 

This is equivalent to the rotations on the unit circle defined by the map 

/cos(27ra;) — sin(27ra;)\ /cos(27r2;)\ 
ysin(27ra;) cos(27ru;) j \sui{2ttx) J ' 

The correspondence between the sequences x„ = M{xn-i) and x„ = i?x„_i is evident after 
writing = (cos(27ra:n), sin(27rxn)). In spite of the rigorous mathematical equivalence between the 
translations on the torus and the rotations on the unit circle, there is not such similitude in the 
numerical evaluation of these maps, as will be explained in the following paragraph. 
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For map Q the divergence of orbits and the reversibihty error have a hnear growth, basically 
due to the the fact that uj* ^ LxJ^ as shown theoretically for the global error in Eq. However 
there are architectures and/or compilers in which the reversibility error may be zero for this simple 
map that only requires the computation of sums and modulus operations. Without knowing the 
precise implementation of these peculiar operations it is impossible to establish a priori with which 
compiler and in which architecture we may obtain this result. We believe this anomaly is due to the 
peculiar arithmetic operations involved and that it does not occur for a generic map. To support 



this claim we have checked that the reversibility error never vanishes for the map ( 18 ) which involves 
multiplications and the evaluation of trigonometric functions. 

In Figure [l]-left we compare the divergence and reversibility error for the torus translation with 
UJ = ^/2 — 1 and xq = 0.7. We can see that both errors approximately satisfy the same expression: 
in log-log scale it is close to a straight line with unitary slope. This implies that both quantities 
are linear functions of time. We have checked that the behavior is similar for almost every initial 
condition and frequency but in some cases the fluctuations around the average linear growth have 
a larger amplitude. Figure [l]-right shows that for the map (18) there is an analogous behavior. 
However, in this case the straight line that corresponds to orbit divergence is a bit shifted upwards 
with respect to the one of reversibility error. This is telling us that the orbit divergence's linear 
growth has a bigger slope. 




100 1000 10000 100000 1e+06 1e+07 



100 1000 10000 100000 1e+06 1e+07 



Figure 1: Comparison between divergence (red) and reversibility (green) for both torus translations 
(left) and rotations (right), using a frequency a; = ^2 — 1 and initial condition xq = 0.7. In all four 
cases the average slope in log-log scale is nearly one. 

The fluctuations Wn have variance < > computed respect to an ensemble of initial conditions 
Xq which first grows linearly, but rapidly saturates to a constant value with small oscillations ( see 
Figure 7 in (Turchetti et al. 2010b[ ) ). In this respect it is quite different with respect to a random 
perturbation e^n where are independent variables, since in this case we have Wn = + - • •+'Cn and 
the variance is given by < w"^ >= n < > for any value of n. This different behavior is reflected 



in the decay of fidelity (Vaienti et al. 2007) which follows a power law for round off errors and 



an exponential law for random perturbations (Turchetti et al. 2010a I , (Turchetti et al. 2010b | 



This shows that the round off errors decorrelate very slowly unlike the random errors which are 
uncor related. 
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For chaotic maps, characterised by an exponential increase in the distance between two nearby 
orbits, the reversibihty error due to round off has the same growth. Bounded two dimensional chaotic 
maps, like the cat map, exhibit this behaviour during a short time before reaching saturation. 



5 The case of the standard map 

As an example of generic 2D map we have chosen the standard map in a torus: 

yn+i = yn + A sin(x„) mod 27r; a;„+i = x„ + mod 27r. (19) 

For very low values of A the divergence of orbits has an average growth linear with n as for the 
translations on the torus and they depend weakly on the initial condition. As A is increased a power 
law with a > 1 is observed and the dependence on the initial conditions becomes appreciable. 
Reversibility error also shows this behavior for A << 1. 

For A approaching one we have coexistence of regular and chaotic regions so that the domain 
is splitted into several islands of stable orbits and a chaotic sea. In figures [2| [3j |4j [5] and [6] we 
show, respectively, the value of the reversibility error (jsl), the orbit divergence ([7|), the finite time 



mLCE (12), the time average of the MEGNO (15) and the SALI (17) obtained iterating eq. (19). 



The standard map with A = 0.971635 has been iterated n = times to compute all the 
dynamical indicators. Even if n < 10'^ might be used to describe the chaotic region, the chosen 
value n = 10^ allows to highlight differences within regular regions where the growth is very slow. 

The pictures in the left side show the value of the dynamical indicators in a chromatic scale 
for a grid of 500 x 500 initial conditions in the two dimensional torus whereas the right pictures 
show the corresponding graph for a fixed value of the action variable (y = 0.3) corresponding to 
the horizontal line in the left figures. Except for mLCE we have used the natural logarithm of the 
absolute value for all the dynamical indicators. 

Figure |2] and Figure |3] were obtained by taking into account the error only in the action variable. 
It is evident that both the reversibility error and the orbits divergence (with respect to the exact 
one) discriminate regular from chaotic orbits. Minor differences exist within the islands of stability: 
the orbit divergence approaches the minimum value close to the centre, whereas the minimum of 
the reversibility error occurs on a line crossing the islands possibly due to a mechanism similar to 



the one reported by Barrio et al. (2009) which analysed spurious errors for variational methods. 

The plots obtained for mLCE and the MEGNO (figures |4] and [S] respectively) show no structure 
within the resonance islands whereas variations and fluctuations appear in the chaotic regions as 
for the previous indicators. 

Figure [6] shows the logarithm of the SALI value. Due to the fact that this indicator converges to 
zero extremely fast for chaotic orbits, we have used a cut off value at SALI = 10~^^. The presence 
of this cut-off erases any structure within the strongly chaotic region. 

To summarise, all the indicators discriminate regular and chaotic regions but their sensitivity in 
these regions is different . 

Another aspect to consider when comparing the efficiency of chaos indicators is the computa- 
tional cost. Each variational method needs to iterate both the map and the tangent map forward 
for n steps. The tangent map is the computationally most expensive since it needs the evaluation 
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of the Jacobian matrix at every step. When computing SALI, two deviations vectors must be si- 
multaneously computed. We also use two deviation vectors in evaluating MEGNO, selecting at 
every step the one that stretches more, in order to reduce the probability of having a vector almost 
orthogonal to the most unstable direction. In the case of the computation of the mLCE we have 



used the orthogonalization algorithm developed by Benettin et al. (1980). On the other hand, we 



notice that the reversibility error method requires only the iteration of the map whereas the orbit 
divergence method requires the iteration of the single and double precision (or double and higher 
precision) map. As a consequence, the computationally most economic method is the one based on 
the reversibility error which does not require any algorithm except the evaluation of the map itself. 
Typically, the relevant information can be extracted from a computation in single precision. 




Figure 2: ln(i?xQ3) displaying the dynamical structure of the phase space of the Standard Map 
(left) and for torus section y = 0.3 (right). 
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Figure 4: mLCE(lO^), displaying the dynamical structure of the phase space of the Standard Map 
(left) and for torus section y = 0.3 (right). 




Figure 5: ln(y(10^)), displaying the dynamical structure of the phase space of the Standard Map 
(left) and for torus section y = 0.3 (right). 
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Figure 6: ln(SALI(10^)) displaying the dynamical structure of the phase space of the Standard 
Map (left) and for torus section y = 0.3 (right). 

From now on, within this section, we will focus on some ensemble quantities in order to compare 
the effect of round off errors with the effect of random perturbations in the standard map. For the 
standard map with A = 10-"^ ah the orbits are regular and we follow the evolution of an ensemble 
of 10001 initial conditions randomly chosen in {x,y) € [1.5, 1.5 + 10~^] x [vr, tt + 10~^]. For each 
iteration we compute the variance of Rn in action (fr^^) and angle (cTa;^) variables (see in Figure]^ 
left) and compare it against the same quantities of a double precision orbit stochastically perturbed 
with a uniform uncorrelated noise of amplitude 10""^ (shown in Figure [Tj-right). In the latter case we 
have found that and ax'^ grow according a power law with exponents one and three respectively 
up to numerical uncertainties. 

It is interesting to compare the two pictures in Figure [7] with each other. For n large we see that 
the behaviors of variances due to round off and random perturbations are similar. The presence of 
a transient for the reversibility case is very likely due to the initial presence of correlation. 

In the fully chaotic regime for the standard map A ^ 1 the random perturbation produces very 
similar results to the round off error as shown in Figure [8] for A = 10 since in presence of a chaotic 
dynamics the round off error correlations are lost. 
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Figure 7: Evolution of the variance of a distribution of reversibility errors in action (red) and angle 
(green) variables for the standard map with A = 10~^. Left: round off error. Right: stochastic 
uncorrelated perturbations of amplitude 10^''. 



Figure 8: Evolution of the variance of a distribution of reversibility errors in action (red) and angle 
(green) variables for the standard map with A = 10. There is evident equivalence between round 
off error (left) and stochastic uncorrelated perturbations of amplitude 10~^ (right). 



The results for the standard map in the small A regime can be compared to the variances for 
the skew map: 

Vn+i = Vn mod 1; x„+i = Xn + Vn+i mod 1. (20) 

to which the standard map in equation (19) reduces for A = 0. If the skew map is randomly 
perturbed: 

Vn+i = yn + e^n mod 1; x^+i = x„ + y„+i + exn modi, (21) 
where ^„ and Xn are random uncorrelated variables the growth of the variances o"^ and cj^. follows a 



linear and cubic law respectively ( Turchetti et al.\ 2006). The random perturbation of a standard 
map with a small value of A (shown in Figure [Tj-right) shows exactly the same growth. 

In the case of the round off error the behavior of the skew map with respect to the standard 
map with a small value of A is quite different. Indeed, the round off error affects the skew map 
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just as the translation on the 1-D Torus: it was observed that the global error grows linearly and 
the variance saturates at a very small value with respect to the size of the torus (see Figure 7 in 
Turchetti et al. (2010a |). For the standard map, the coupling between action and angle, even for 



very small A, causes a growth of the variance of the fluctuations due to round off as shown in Figure 
[^left. As a consequence the effect of round off in a very weakly perturbed standard map is similar 
to a random perturbation. 



6 A 4D Map 

In this section we show how either the reversibility error or the divergence of orbits can be used 
to analyze the resonance structure of four dimensional non integrable maps. As an example we 



consider a symplectic nearly integrable map extensively used in the literature (see Guzzo & Lega 



(2004|). This map is defined as: 



6n+l 


— 


+ In 




= 4>n 


+ Jn 


In+1 


= In 


dV{dn+l,(l)n+l) 
^ dOn+l 


Jn+1 


— Jn 


dViOn+lAn+l) 



i^n+l 



(22) 



where V = l/(cos(0) + cos((/>) + 2 + c), with c > and fi is the perturbation parameter. We 
have taken a grid of 1146 x 1146 initial conditions with actions (/, J) € [0, 3.6] x [0, 3.6] and a 
fixed pair of angles, namely (9,(1)) = (0.5, 0.5) and computed Rn, A„, mLCE(n), Y(n) and SALI(n) 
for n = 10^. Associating these values to each initial condition while using a chromatic scale we 



have performed figures [9j [TO] and 11, where the parameter values c = 2 and fi = 0.6 have been 
used. As we did in the previous section, in all except for the mLCE we have used the natural 
logarithm of the absolute value of the concomitant indicator and again we have used a cut off value 
at SALI = 10^^®. Figure [ojwas done taking into account the Euclidean error in only the action 
plane. As it happened for the standard map both reversibility error and orbit divergence present 
the same order of magnitude. The resonance web appears in every figure and its structure is the 
same. From a qualitative viewpoint no substantial differences are found and one might conclude 
that the dynamical information extracted from the reversibility error is the same as for the other 
four dynamical indicators we have considered. In the case of mLCE, the highest values are located 
on the diagonal I = J even if they are not well visible in the plot. This is anyway coherent with 
the fact that on the diagonal we find the highest instability. 

The frequency map analysis was not presented for comparison because it is computationally heavier. 
The error reaches its highest values in small neighbourhoods of single resonance lines, because 
it detects the perturbed separatrices of the pendulum models that one could associate with the 
resonances, and in relatively large neighbourhoods of the intersections of resonances, where the 
dynamics are widely non integrable, as shown by the computation of the interpolating Hamiltonian. 
These results are certainly not exhaustive but show that the behavior of the reversibility error is 
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strictly related to the divergence of orbits and consequently it is very weak in the integrable regions 
where it does not have a diffusive character as for a random error. 





7 Conclusions 

We have examined the orbit divergence and the reversibihty error in order to determine the effect 
of round off error for invertible maps. The knowledge of the exact map is not required to compute 
the reversibility error and the results obtained are about the same with respect to the case in which 
we compute the orbit divergence. 

By choosing an ensemble of initial data we have examined the statistics of the fluctuations due 
to round off with respect to the time average of the error. There are different behaviours according 
to the degree of chaoticity that characterises the ensemble. For chaotic orbits the variances have an 
exponential growth similar to the one observed for random perturbations and correspondingly the 
decay of fldehty is super-exponential in both cases. For regular or quasi-regular orbits differences are 
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Figure 11: SALI(IO^) for map (22) using c = 2 and ^ = 0.6. 



observed between the effects generated by round off errors and random perturbations. For a quasi- 
integrable map the random perturbations produce a growth of variances hnear and cubic for actions 
and angles, respectively. On the other hand, the round off errors produce an initial transient that 
possibly lasts the time that the round off errors of the ensemble orbits need to decorrelate. After 
this initial transient, the behaviour of variances affected by round off and noise are similar to each 
other. Finally, for a regular map such as the translation of the ID torus, the round off variance, after 
a linear growth, saturates before the distribution of errors fills all the torus, unlike in the randomly 
perturbed case where the growth is always linear until saturation of the full torus. The fidelity has 
a power law decay for round off whereas it decays exponentially for random perturbations. This is 
a clear signature of the correlation between the errors due to round off. 

To conclude, the reversibility error provides basically the same information as the divergence 
of orbits and it is easily accessible from a computational view point. This is due to the facts 
that it does not require the solution of the variational equations and that both the forward and 
backward orbits can be computed using single precision. When the initial point is varied and the 
number of iterations is kept fixed the reversibility error due to round off provides an insight on the 
dynamical structure of the map. For the standard map the reversibility error provides a picture of 
the dynamical behavior of the map where not only large scale features but also small scale details 
can be detected. In the case of a 4D symplectic map the reversibility error in action space provides 
a similar picture where the resonance web and the nearby regions of weakly chaotic motions can be 
easily highlighted. 

Even if no really new information with respect to the standard indicators is provided by the 
reversibility error, we point out that this type of analysis takes into account not only the dynamics of 
the map but also the unavoidable effect of finite accuracy due to round off. Moreover, by increasing 
the accuracy of numerical computations, the time on which the reversibility test is computed can 
be increased and finer details on the phase space structure can be detected. 
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